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ABSTRACT 

The large-scale structure of high-redshift galaxies produces correlated anisotropy in the far-infrared 
background (FIRB). In regions of the sky where the thermal emission from Galactic dust is well below 
average, these high-redshift correlations may be the most significant source of angular fluctuation power 
over a wide range of angular scales, from ^7' to ^--^3°, and frequencies, from ^^400 to ^1000 GHz. The 
strength of this signal should allow detailed studies of the statistics of the FIRB fluctuations, including 
the shape of the angular power spectrum at a given frequency and the degree of coherence between FIRB 
maps at different frequencies. The FIRB correlations depend upon and hence constrain the redshift- 
dependent spectral energy distributions, number counts, and clustering bias of the galaxies and active 
nuclei that contribute to the background. We quantify the accuracy to which Planck and a newly 
proposed balloon-borne mission EDGE could constrain models of the high-redshift universe through the 
measurement of FIRB fluctuations. We conclude that the average bias of high-redshift galaxies could be 
measured to an accuracy of ^1% or, for example, separated into 4 redshift bins with ~10% accuracy. 

Subject headings: cosmology: theory - cosmology: observation - cosmology: far infrared background - 
cosmic microwave background - galaxies: formation - galaxies: evolution 



1. INTRODUCTION 

The discovery of the cosmic far-infrared background 
(FIRB; Puget et al. 1996; Fixsen et al. 1998; Dwek et 
al. 1998; Schlegel et al. 1998; Lagache et al. 1999) and 
the determination of its spectrum have resulted in a new 
probe of structure formation in the high-redshift universe 
(e.g. Guidcrdoni et al. 1998; Blain et al. 1999a; Haiman & 
Knox 2000). A compelling explanation for the background 
is that it results from the thermal emission of interstellar 
dust associated with high-redshift galaxies and heated by 
the internal optical and ultraviolet (UV) radiation from 
stars and to a lesser extent active galactic nuclei (AGN) 
(Stecker et al. 1977; Bond et al. 1986). Thus, the basic 
properties of the FIRB are sensitive to the role of galaxy 
formation and subsequent evolution. 

Observing the FIRB at higher angular resolution and 
higher sensitivity than was done by the Far Infrared Ab- 
solute Spectrometer (FIRAS, Mather et al. (1999)) instru- 
ment on the Cosmic Background Explorer {COBE) satel- 
lite should reveal the presence of correlated fluctuations 
resulting from the correlations in the galaxies contribut- 
ing to the background. Haiman and Knox (2000, hereafter 
HKOO) used simplified semi-analytic models for the origin 
of the background flux to show that these correlated an- 
isotropics have an amplitude that is roughly 10% of the 
mean, a level detectable with current technologies. The 
purpose of the present paper is to (1) spell out and quan- 
tify what one could learn from detailed observations of 
the FIRB correlations, (2) show that such detailed studies 
are possible due to the absence of other strong sources of 
fluctuation power in the relevant range of frequencies and 
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angular scales, and (3) show that FIRB anisotropy ob- 
servations provide a powerful complement to direct high- 
angular resolution observations of the individual sources. 

As we discuss in detail later, the background at different 
frequencies is composed of sources from differing ranges of 
redshifts. FIRB maps at multiple frequencies are there- 
fore not expected to correlate perfectly with each other. 
We show how the shape of the FIRB angular power spec- 
trum at different frequencies and the correlations between 
the maps can be used to determine physical properties of 
contributing high-redshift sources, in particular the prod- 
uct of the infrared emissivity of the sources and their bias 
relative to the dark matter density field. With assump- 
tions about the source biasing and spectra, the infrared 
emissivity as a function of redshift can be converted to 
a measurement of the energy-production history of the 
universe, complementing current approaches involving op- 
tical and UV observations (e.g. Madau 1997; Madau et 
al. 1998). Conversely, one could combine FIRB anisot- 
ropy measurements with other measures of the emissivity 
density — e.g, from deep sub-millimeter surveys followed up 
by redshift determinations — to determine the bias of the 
sources. 

Other important sources of fluctuation power over the 
relevant range of frequencies and angular scales are the 
cosmic microwave background (CMB), the shot noise due 
to the discrete nature of the FIRB sources, and thermal 
emission from dust in our own Galaxy. The CMB prevents 
measurements of the FIRB at frequencies less than about 
200 GHz. Even at 200 GHz the angular power spectrum 
of the CMB is over an order of magnitude larger than the 
FIRB's, but with a sufficiently sensitive lower frequency 
CMB map one could use the well-known spectral depen- 



dence of the CMB to subtract it with high precision. We 
expect the correlated FIRB fluctuations to dominate the 
shot-noise contribution at multipole moments t < 1500, 
or angular scales > 7'. Thermal emission from dust in our 
own galaxy could easily be confused with FIRB anisot- 
ropy, but, as we show below, the dust fluctuation power 
has a strong spatial dependence. There are regions of suf- 
ficiently low dust fluctuation power that the FIRB fluc- 
tuations are dominant at ^ > 60 (angular scales < 3°) at 
low frequencies. This critical multipole moment slowly in- 
creases as the frequency increases towards v ^ 1000 GHz 
and then rapidly increases beyond there, as the Milky Way 
dust appears hotter than the high-z dust contributing to 
the FIRB. Between these frequency and multipole moment 
limits, the FIRB correlations can be studied without the 
need for aggressive foreground subtraction. 

Our work has been largely motivated by the proposed 
balloon-borne Explorer of Diffuse Galactic Emissions^ 
(EDGE), which will survey the sky in 10 frequency bands 
from 150 GHz to 1290 GHz, with angular resolution rang- 
ing from 14' to 6' full-width half-maximum (FWHM). 
Later in the decade, the Planck surveyor"^ should provide 
high-quality maps of correlated FIRB anisotropy at 217, 
320, 545 and 850 GHz. We forecast the errors on the 
scale- and redshift-dependent product of bias and emissiv- 
ity that could be reconstructed from EDGE and Planck ob- 
servations. We also discuss how these quantities could be 
disentangled from uncertainties in the cosmological model 
and the emission spectrum of the sources. 

Both sets of observations can reconstruct the product 
of emissivity and bias to ^ 10% in 4 redshift bins with 
the highest one extending from 2; = 2 to z = 4. Indeed, 
the constraints are ~ 1% if one considers smoother models 
for the bias-weighted emissivity. This level of precision 
is hard to duplicate by directly observing the clustering 
properties of the sources. It may be possible to resolve the 
background as a function of redshift in several frequency 
bands in the next decade using, e.g., BOLOCAM (Glenn 
et al. 1998) on the Large Millimeter Telescope (LMT) or 
the Atacama Large Millimeter Array** (ALMA) and doing 
follow-up optical redshift determinations. However, we 
show that determining the clustering of the sources to < 
10% accuracy at large enough scales to be in the linear 
regime requires many thousands of sources with follow-up 
redshifts of very faint galaxies. 

Other missions may also detect the correlated FIRB an- 
isotropy. The TopHat experiment^, with its 660 GHz chan- 
nel, may provide an early opportunity for studying the 
FIRB anisotropy. The balloon-borne Bolometric Large 
Aperture Sub-millimeter Telescope® (BLAST), designed 
to detect sources at the bright and rare end of the source 
count distribution, may also be able to measure fluctua- 
tions in the diffuse background. 

Indeed, detectable FIRB correlations may even be lurk- 
ing in existing data sets. Using the ISOPHOT instru- 
ment on ISO^, the shot noise of the FIRB has been seen 

■^ http ; / /topweb . gsfc . nasa. gov 

3http;//astro.estec.esa.nl/Planck/; also, ESA D/SCI(6)3. 

^http;//www. alma.nrao.edu 

''littp://topwcb. gsfc. nasa. gov/ 

®littp://www. hep. upcnn.edu/blast/ 

^sce http://isowww.cstcc.csa.nl/ for details 



at 170y^m (1760 GHz) but the galactic dust emission, at 
least in the field observed, has obscured the FIRB cor- 
relations (Lagache & Puget 2000). The correlations may 
be detected upon further analysis of the ISO data (after 
some dust cleaning), or in the BOOMERanG 420 GHz 
data (de Bernardis et al. 1999), or even possibly in the 
COBE/FIRAS data. Kashhnsky & Odenwald (2000) claim 
a tentative detection of correlated fluctuations in the near- 
infrared background using COBE/DIRBF, data shortward 
of 10 yum, with somewhat greater amplitude than that pre- 
dicted by Jimenez & Kashlinsky (1997). 

The remainder of our paper is organized as follows. In 
§ 2, we review what is known about the sources of the 
FIRB and about the mean spectrum. In § 3, we describe 
both the correlated and shot-noise contributions to the 
anisotropy. For the correlated component we demonstrate 
the dependence of the angular power spectrum on the fluc- 
tuation power in the FIR emission at different redshifts. 
In § 4, we describe the fluctuation power in the sky at the 
frequencies of interest from the other relevant components: 
thermal emission from dust in our own Galaxy, the CMB, 
and the FIRB. In particular, we emphasize the wide dis- 
tribution in amplitudes of dust power spectrum across the 
sky. In § 5, we discuss qualitatively the reconstruction of 
various interesting physical quantities from FIRB obser- 
vations, and then forecast EDGE and Planck errors on a 
reconstruction of a scale- and redshift-dependent product 
of bias times emissivity. We discuss our results in § 6, with 
an emphasis on the complementary nature of point source 
observations, and summarize our conclusions in § 7. 

Although we maintain generality in all derivations, we 
illustrate our results with the currently favored AGDM 
cosmological model. Our choices for cosmological parame- 
ters are n^ = 0.30, f^;, = 0.05, f^A = 0.65, h = 0.65, n = 1, 
and GOBE normahzation of 6h — 4.2 x 10^^ (Bunn & 
White 1997). We use the fitting formula for the transfer 
function given by Eisenstein & Hu (1999). This model has 
mass fluctuations on the 8h~^ Mpc scale in accord with the 
abundance of galaxy clusters as = 0.86 (Viana & Liddle 
1999). 

2. FIRB MEAN 

In this section, we review what is known about the 
FIRB mean, both from the FIRAS determination of its 
spectrum and from observations of sources — most notably 
at 850/zm by the Sub-millimetre Common User Bolome- 
ter Array (SCUBA) camera (Holland et al. 1998) on the 
James Clerk Maxwell Telescope (JCMT). We emphasize 
the connection, created by the redshifting of the peak of 
FIR emission, between studying the background spectrum 
at low frequencies and determining the FIR emissivity at 
high redshifts. 

2.1. The Spectrum 

A grey body with power-law emissivity has been fit to 
the spectrum of the mean FIRB as determined with the 
COBE/FIRAS data (Fixsen et al. 1998, hereafter F98): 

I,^To(iy/i^o)''B,(T) (1) 

where B^(T) is the Planck function and tq is the optical 
depth at t/Q = 3 THz. F98 find tq = (1.3 ± 0.4) x 10"^ 
T — 18.5 ± 1.2K and emissivity power-law index a = 
0.64 ± 0.12; this spectrum is plotted in Figure 1. 
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Fig. 1. — The spectrum of the FIRB predicted in the HKOO stan- 
dard model (solid curve). The dashed curves show the measure- 
ments with ±1(7 uncertainties from F98. The light dotted curves 
are grey-body spectra with a = 2 and T = 8 K (left curve) and 
T = 14 K (right curve). Their optical depths and temperatures 
are chosen to match the low frequency and high frequency sides of 
the HKOO spectrum. The HKOO spectrum is a sum of many such 
grey-body spectra. 



HKOO considered a variety of models to describe the 
FIRB spectrum. For the present calculations, we use only 
the "standard model" from HKOO, which we refer to sim- 
ply as the HKOO model. Its spectrum is shown in Fig. 1, 
with the amplitude reduced by 40% relative to HKOO to be 
in better agreement with the FIRAS data. We apply the 
reduction throughout this paper, thereby reducing the an- 
gular power spectra by a factor of 2. The HKOO model for 
the FIRB mean assumes that the UV emissivity at each 
rcdshift is proportional to the star-formation rate (SFR) 
as estimated by Madau (1999), that dust production is 
proportional to the SFR, and that the dust has the opti- 
cal properties of the Draine and Lee (1984) model. With 
the optical properties and radiation backgrounds specified, 
the dust temperature follows. The two proportionality 
constants for the UV emissivity and dust production are 
set to produce a spectrum that is consistent with the F98 
determination. 

On the Rayleigh-Jeans side of the spectrum, Draine and 
Lee (1984) dust is a grey body with emissivity index a — 2. 
Thus the HKOO spectrum is a sum of such grey bodies with 
varying apparent temperatures and optical depths, where 
apparent temperature is defined as Tapp = T/{l + z) where 
T is the physical temperature. The HKOO spectrum can- 
not be fit by a single grey body with a = 2. We demon- 
strate this in Fig. 1 with the a — 2 grey bodies plotted 
as dotted curves. For the right-most dotted curve, we 
chose the temperature and optical depth to fit the curve 
at high frequencies. Moving towards lower frequencies. 



this Tapp = 14 K grey body quickly falls below the HKOO 
spectrum. To fit the spectrum, one must add in colder 
components, such as the Tkpp = 8 K one shown. Clearly, 
in this particular model, there is not much need for even 
colder components. Since the HKOO model has dust with 
r « 24 K in the rcdshift range 1 to 4, this indicates that 
most of the FIRB in this model comes from z < 2. Note 
that the F98 fit to the FIRAS data is even softer than the 



HKOO model, possibly due to the presence of even colder 
(higher rcdshift) components, although it should be kept 
in mind that the uncertainties in this slope at low frequen- 
cies are quite large. 

In principle, one can take the spectrum of the back- 
ground and, assuming the emissivity index, reconstruct 
the optical depth of each temperature component. This is 
the essence of the idea recently discussed and implemented 
by Gispert et al. (2000). By assuming the shape of the av- 
erage SED, one can recover the FIR luminosity density as 
a function of rcdshift, which can be related to the SFR 
using certain assumptions. 

Due to the large size of the FIRAS uncertainties, es- 
pecially at frequencies less than ^ 300 GHz, this recon- 
struction can not be done accurately. Furthermore, there 
are no proposed missions for improving upon FIRAS 's de- 
termination of the FIRB spectrum, which wc understand 
to be a very difhcult task. Fortunately, one can make 
quite similar use of a measurement of the spectrum of the 
FIRB anisotropy — although with the added complication 
that one is not reconstructing the luminosity density, but 
the luminosity density weighted by the clustering bias of 
the sources. A differential measurement of the spectrum 
should be able to measure to lower frequencies due to the 
fact that the amplitude of the CMB anisotropy is about 
10^'"' of the CMB monopole. This is especially interesting 
since the limiting rcdshift (coldest component) one can 
reach is determined by the lower frequency limit of the 
spectrum determination. We turn to a description of the 
FIRB anisotropy in § 3. 

2.2. Resolving the Background 

The mean brightness of a background due to point 
sources depends on the number density of sources as a 
function of limiting flux, N{>S): 



{S'^dN/dS)d\nS. 



(2) 



Notice that the bulk of the background is contributed by 
sources near the maximum of S'^dN/dS. We therefore 
show this function as predicted by the semi-analytic mod- 
eling of Guiderdoni et al. (1998, hereafter G98) at several 
frequencies in Fig. 2. G98 calculated predictions for a 
range of models; we use their model E predictions, which 
are the ones in best agreement with the observations avail- 
able at 850/im. 

Considerable progress has been made towards resolv- 
ing the FIRB into discrete point sources (Hughes et al. 
1998; Barger et al. 1998; Eales et al. 1999; Small et al. 
1997; Holland et al. 1998; Blain et al. 1999b; Barger et 
al. 1999a; Puget et al. 1999; Barger et al. 1999b; Blain 
et al. 1999c). In particular, the SCUBA camera (Holland 
et al. 1998) on the JCMT has been used to identify point 
sources at 850/im that account for a large fraction, ~ 50%, 
of the FIRB at these wavelengths (Barger et al. 1999a). As 
shown in Figure 2, the SCUBA detections are in the range 
of ~ 1 mjy to ~ 10 mjy, presumably near the peak of 
S'^dN/dS. The result of fitting a double power-law form 
for iV(> S) to the SCUBA data, with the constraint that 
it reproduce the FIRAS-determined mean (Barger et al. 
1999a), results in the light solid line in Figure 2. 

FoUowup observations of SCUBA locations at radio, 
optical and infrared wavelengths have suggested possible 
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Fig. 2. — The contribution to the FIRB mean from each log- 
arithmic interval in S, according to model E (G98) at 300 GHz, 
353 GHz and 400 GHz. The light dashed line is a fit to the SCUBA 
data, which covers the ~ 1 mJy to ~ 10 mJy range (Barger et al. 
1999a). Future observations near these frequency ranges will go 
deeper and broader. The left— hand limits for BOLOCAM, observ- 
ing at 270 GHz, are the 5(7 sensitivities after a 10-hour pointing 
at a single field-of-view. Both of these pointing should result in 
the detection of about 60 sources according to model E (G98). The 
left— hand limit for BLAST is the 5(7 limit for the deepest of their 
planned surveys, which takes 2 days of a 10— day long— duration bal- 
loon flight. The vertical location of the indicated survey ranges is 
arbitrary. 



counterparts, albeit with a relatively low identification 
rate. Typical sources are identified as galaxies with mod- 
erate to massive star formation rates. Optical redshifts, 
the sub-millimeter to radio spectral index (e.g. Carilli & 
Yun 1999) and the detection of CO molecular lines (Frayer 
et al. 1999) suggest that most of these sources lie at red- 
shifts ranging from 1 to 4. For a recent review of SCUBA 
results see Small et al. (2000). Continued foUowup obser- 
vations of SCUBA surveys and those that were carried out 
with ISO, such as the European Large Area ISO Survey 
(Efstathiou et al. 2000; Serjeant et al. 2000; Oliver et al. 
2000), will provide important details on the physical prop- 
erties of far-infrared (FIR) emission and the redshifts of 
contributing sources. 

Further progress towards resolving the background in 
wave bands near 850/im (353 GHz) will come soon from 
a bolometer array sensitive to 1100/im (273 GHz) called 
BOLOCAM, now at the Caltech Submm Observatory 
(CSO). When deployed on the LMT, it wiU have suffi- 
cient resolution to resolve considerably more of the back- 
ground. BLAST wiW be able to detect the brightest sources 
at 750/im (400 GHz) and also has detectors at 300 and 
200/xm, to permit the determination of color redshifts. 
The California Millimeter Array (CARMA), an interfer- 
ometer composed of existing BIMA and Owens Valley tele- 
scopes, is expected to carry out observations at 235 and 
345 GHz with ^ arc second resolution that will allow cross- 
identification at optical wavelengths, allowing further fol- 
lowup observations (Carlstrom 2000). Finally, ALMA will 
have high sensitivity and very high angular resolution. In 
a single 10 hour pointing, it will achieve a 5cr noise level of 
0.03 mJy and still be far from the confusion limit. How- 
ever, many such pointings will be required to get large 
numbers of sources, as the field-of-view is only 0.07 square 



arcminutes, over which model E predicts about six sources 
brighter than 0.03 mJy. 

At shorter far-infrared wavelengths, a smaller fraction of 
the background has been resolved. For example, the FIR- 
BACK survey at 170 /xm with ISOPHOT on ISO (Puget et 
al. 1999) has resolved about 10% of the background (Puget 
et al. 1999). Much progress will come in the near future 
from MIPS on SIRTF (160, 70 and 24/xm), SPIRE and 
PACS on FIRST (90 to 500^nJi), and HAWC on SOFIA 
(60, 110 and 200 fim). Information on most of these sub- 
millimeter and far-infrared instruments is summarized in 
Blain (1999d). 

3. FIRB ANISOTROPY 

The FIRB is anisotropic at small scales due to the dis- 
crete nature of the sources, and at larger scales due to 
correlations between these sources. We first discuss the 
correlated anisotropy. 

3.1. FIRB correlations 

The antenna temperature of the FIRB at a given fre- 
quency v and in a given direction n can be written as a 
line of sight integral of the product between the mean FIR 
emissivity and its fluctuation: 



dz—a{z)j{v,z) 
dz 



1 



5j{r{z)n,v,z) 



(3) 



]{v,z) 

where r is the coordinate distance (or conformal time) 
from our location at the coordinate origin and j(z^, z) is the 
mean emissivity per comoving unit volume at frequency v 
as a function of redshift z. In a flat universe, which we 
will assume, r is simply the proper motion distance. We 
will always use comoving wavenumbers and spatial scales. 
We next wish to derive the angular power spectrum C^" 
of the FIRB, which is the Legendre transform of the two- 
point correlation function 

C(n, rii) = (T(fi, i/)T(m, v')) 

'^^^^^"•''Pi{n-rh). (4) 
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This is most easily done by first decomposing the temper- 
ature maps into spherical harmonic multipole moments. 



ai7n{v) 



dfiT(fi,j/)y,™*(n). 



(5) 



With an isotropic random field, the second moments of 
these multipole coefficients arc 

K™,Ma^.m.(^')) = €l.C^^.C'^f■ (6) 

Using equation (3), we can rewrite equation (5 as 



a™{^)^ j dnYr{n)T{n,v) 



where 



|^<5(k)/™B(fc^^)^„»(j^), (7) 

/™(fc^^)= I dz'^W™^{k,i.,z)j,{kr) 

dr 



VK^^^^"(fc,i/,z) / dz—ji{kr) 
J dz 
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and 

W^^^^ (fc, ly, z) = a{z)G{z)j{iy, z)b{k, v, z) . (9) 

Here, G{z) is the linear theory growth function (Peebles 
1980), J£ is the usual Bessel function, and the approxi- 
mation used follows the well-known Limber form (Limber 
1953). The ratio of Gamma functions goes to ^2/1 for 
€>1. 

In equation (7), (5(k) is the linear theory dark matter 
density fluctuation today, 5 = 5p/p. It enters because we 
have assumed, as in HKOO, that the fluctuations in the 
emissivity, 5j / j, are a biased tracer of those in the mass. 
In general, this bias can depend on scale, frequency, and 
redshift, so that in Fourier space we write: 

<5j(k, t^,z) 



6(fc, V, z)(5(k, z). 



(10) 



We ignore the possibility of stochastic bias (Dekel & Lahav 
1999) until the discussion. 

With these assumptions about the bias we can finally 
write the angular power spectrum of FIRB anisotropy us- 
ing the three-dimensional, linear-theory power spectrum 
of dark matter density fluctuations today, PM{k), 

{6(k, z)S*{k', z)) - (2^)3j^(k - k')PM(fc)G2(z) , (11) 

as a line of sight projection 

-iT—a^{z)jb{k, V, z)jb{k, v\ z)PM{k)\k=LG' 



dz 

We have used the Limber approximation which sets k = 
l/r. At £ = 30 the Limber approximation is valid here 
to within 10% of the exact (first order) calculation and 
rapidly converges to the exact value as i increases. Since 
the bias factor enters into the expression for C^''^ multi- 
plied by the mean emissivity, we define the "bias-weighted 
emissivity" jt as 



jh(fc, ly, z) = b{k, ly, z)j{iy, z) . 



(13) 



We do not use the above equations in their full gen- 
erality. We follow HKOO in assuming that the bias is a 
constant, i.e., independent of frequency, redshift and scale. 
We take this bias to be 3, which will roughly match the ob- 
servations of the clustering of z = 3 Lyman-break galaxies 
(Steidel et al. 1998; Giavalisco et al. 1998; Adelberger et 
al. 1998). Bias is likely to increase with redshift (Baugh 
et al. 1999; Blanton et al. 2000), and hence this approx- 
imation may down-weigh the relative importance of fluc- 
tuation power at high redshift. The bias may also be fre- 
quency dependent, which could result from two different 
populations of sources, each with different spectral energy 
distributions and clustering properties. We briefly discuss 
such a possibility below but neglect it for now in order 
to illustrate the nature of the observables in the simplest 
possible case. 

We wish to understand how measurements of Cj^'^ can 
be used to determine the bias-weighted emissivity jb as a 
function of redshift. To that end we display contributions 
to the angular power spectra from various redshift ranges 
in Figure 3. As stated above, we have assumed that the 
bias is redshift-independent and that the mean emissivity 
is that of the HKOO standard model. 

Figure 3 shows that the relative contributions from the 
various redshift ranges depends on the observed frequency 
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Fig. 3. — dCY^ /dz for a power- law 3D power spectrum (solid 
lines) at, from right to left, u = 210, 353 and 1090 GHz, all nor- 
malized to unity at their maximum. For any power-law 3D power 
spectrum, these are independent of I. The three curves in the bot- 
tom panel show the same quantity, using a ACDM power spectrum 
at 353 GHz, for i = 30, 100, and 500 respectively 



and on the shape of the matter power spectrum. In 
the top panel, we consider a power-law power spectrum, 
P{k) ex fc^^, and show the shape of the redshift break- 
down dCY^/dz for three frequencies, ly = 210,353 and 
1090 GHz. These are the integrands of equation ( 12). 
Note that lower frequencies probe to higher redshifts, as 
expected. For any power-law power spectrum, the shape 
of dCi'^ /dz as a function of z is the same for all values of 
l. 

The non-power-law nature of physical power spectra 
lead to different shapes of dCi" /dz for different values of 
i. In the lower panel of Figure 3, we plot three different 
i at fixed frequency for the ACDM model. Higher £ map 
to smaller scales, where the spectral index of the power 
spectrum is more negative. This means that higher red- 
shifts, where the spatial scale corresponding to a given i. is 
larger, contribute more to the anisotropics. Hence, higher 
d probe higher redshifts. 

Figure 4 shows the contributions to Ci'^ for different 
redshift ranges and frequencies. In each panel, one can 
see the shape dependence of the contribution from each 
redshift interval, with lower redshifts having more low-£ 
fluctuation power compared to the high-redshift intervals. 
This shape dependence aids in the reconstruction of the 
redshift dependence of the fluctuation power. Comparing 
the panels, one should note the more rapid drop-off to- 
wards higher redshift at 850 GHz than at 210 GHz. At 
the higher frequency, there is very little contribution from 
z> 1.5. 

We now turn to the cross-frequency correlations. In Fig- 
ure 5, we show the expected cross-correlations between 
measurements at 390 GHz and those at four other fre- 
quencies (210, 852, 1060, and 1290 GHz). The lowest fre- 
quencies are affected by all the redshifts that have any 
signiflcant fluctuation power; therefore, these frequencies 
weight all of the fluctuations in the same way and are 
therefore highly correlated. Higher frequencies are insen- 
sitive to high-redshift fluctuations because they fall on the 
Wien side of the redshifted spectrum. With only a por- 
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— Contributions to C;"" from several redsliift ranges 
(labeled in lower panel) for 1/ = 210 GHz (upper panel) and 
u = 850 GHz (lower panel). The heavy solid line is the total power 
spectrum observed at z = 0. 



tion of the fluctuations included, the degree of correlation 
degrades. The break frequency, at which this correlation 
departs from unity, is sensitive to the onset of significant 
fluctuation power in FIR galaxies. Because higher £'s are 
preferentially more sensitive to higher z, the diminishing 
importance of high-z at higher frequencies is more pro- 
nounced at smaller angular scale. 

It is important to keep in mind that the redshift- 
dependences of the power spectra, illustrated in figures 
3 and 4, and the cross-correlations shown in figure 5 are 
those of a particular model — the HKOO standard model. 
For this model the emissivity tracks the star formation 
rate as determined by Madau et al. (1998), which is quite 
uncertain. It is possible that the SFR is much higher at 
high-redshift. If one assumes efficient conversion of gas 
into stars in all collapsed halos down to some small mass 
(with virial temperatures of abut 10'' K) then there is a 
peak in the SFR at z ~ 7 (e.g. Haiman et al. 2000) which 
is seen in simulations (Gncdin & Ostriker 1996). Increased 
high-z SFR would lead to a greater fraction of the fluctua- 
tion power coming from high redshifts, and greater decor- 
relation between the high and low frequencies. 

We use the linear-theory evolution of the matter power 
spectrum in our calculations. At z ss 1, the non- 
linear corrections become important at wavenumbers k > 
0.4ft. Mpc^^, which projects to i > 800. At lower redshifts, 
this non-linear scale moves to smaller £, but the FIRB con- 
tributions drop as well (see the higher i curve in Fig. 3). 
At ^ « 1000, non-linear corrections produce a 20% in- 
crease in the FIRB angular power spectrum. Because this 
is a small difference that can be accurately included for 
any chosen cosmology, neglecting the non-linearity won't 
affect our quantitative results. However, we do restrict our 
analysis to i? < 1600 both to avoid the deeply non-linear 
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Fig. 5.— Cross-correlations C^"' /Jcj^^'Cj^'''' vs. l between the 

u = 390 GHz channel and four other frequencies u' = 210, 852, 
1060, and 1290 GHz. 



regime and to eliminate the need for detailed shot noise 
subtraction. 

3.2. Shot Noise 

At smaller angular scales, the fact that the FIRB comes 
from individual, well-separated sources introduces anisot- 
ropy simply due to the Poisson sampling of the density 
field. This so-called shot noise can be calculated from the 
source counts by 



yOShot 



{S^dN/dS)dlTiS. 



(14) 



Note that this is more heavily weighted towards the bright 
end of the distribution than the mean brightness was 

For a Euclidean distribution (appropriate for galaxies 
at low redshift), dN/dS ex S^"^'^ and so the integrand di- 
verges at the bright end. However, in the case of the FIRB 
(at least as modeled by G98 model E), the enormous num- 
ber of faint, high-redshift sources cause the source counts 
to become steeper than S~'^^^ at fainter flux levels. Even if 
we remove sources only down to the very conservative flux 
level S'cut such that N{> Sent) ~ 1 per 47r steradians, the 
FIRB shot noise is still dominated by the fainter, high- 
redshift population. Therefore, we take our shot noise es- 
timates from the G98 model E counts, assuming only a 
cut at one source per sky and avoid any complications of 
source removal. The number of sources detected at greater 
than 5(7 by EDGE, according to model E, are given in Ta- 
ble 1 . The removal of these sources from the map does not 
significantly reduce the shot noise. 

4. FOREGROUNDS AND BACKGROUNDS 

There are celestial sources of emission in the far infrared 
and sub-millimeter other than the FIRB. To detect the 
correlations of the high-redshift galaxies that make up the 
FIRB, wc need to be able to separate this signal from all 
the other sources. The two most important contributions 
come from the CMB anisotropics and thermal emission 
from dust in our own Galaxy. Radio foregrounds such 
as synchrotron, Bremsstrahlung, and rotational emission 
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of the cleanest regions, respectively; and the thick (square) dotted curves with spacings of SI = 50 are errors expected on the determination 
of the FIRB power spectrum from a 10-day flight of EDGE. 



from dust do not appear to be significant at z^ > 200 GHz. 
Reviews of wfiat is known about galactic dust emission and 
other CMB foregrounds can be found in dc Oliveira-Costa 
& Tegmark (1999). 

Figure 6 shows the predicted power spectra of the com- 
ponents at 4 different EDGE frequencies. The FIRB an- 
gular power spectrum is shown by the thick dashed curves 
for the standard, constant bias model of HKOO. Unlike 
in all previous figures in this paper, these curves include 
a shot-noise contribution, computed from the G98 model 
E, as well as the effects of the non-linear evolution of the 
matter power spectrum. Other lines show the contribu- 
tions of the CMB and Galactic dust, which we will now 
discuss in detail. 

4.1. CMB 

The thick solid curves in Figure 6 show the CMB power 
spectrum in our flat ACDM model, consistent with the 
COBE/Y)MR, Boomerang and Maxima measurements (de 
Bernardis et al. 2000; Hanany et al. 2000; Jaffe et al. 2000). 
As one can see from this figure, the CMB anisotropics are 
the dominant source of fluctuation power on the sky over 
a substantial range of £ at low frequencies. The FIRB 
anisotropics would be swamped by such a signal. For- 
tunately, the spectrum of the CMB anisotropics is well- 
known, so that one can use observations at yet lower fre- 
quencies (where the contrast between the CMB and the 
FIRB is even starker) to clean the CMB signal out of the 
FIRB maps. 

To demonstrate the ability of a particular experiment 
to subtract the CMB, we have forecasted errors on the 
determination of the FIRB power spectrum for 4 of the 
EDGE channels. The results, in bins of 61 — 50, are plot- 
ted in Fig. 6. For each of the error forecasts at 210, 320, 
590 and 1060 GHz, we use only two channels — the channel 



in question plus the 150 GHz channel (as a CMB moni- 
tor). Note that the 590 and 1060 GHz channels require 
no CMB cleaning, and therefore the FIRB measurement 
is sample-variance limited until the beam becomes impor- 
tant at i> 1000. At 210 and 320 GHz, the errors increase 
above sample variance at much lower i, due to imprecise 
cleaning of the CMB. Even at these low frequencies, the 
CMB can be subtracted sufficiently accurately to allow for 
SCi/Ci < 1 in bins of SI = 50 up to ^ - 1000. The more 
sensitive and higher angular resolution 100 GHz channel 
on Planck allows for even better CMB subtraction. More 
details on these power spectrum error forecasts are given 
in section 5.2. For now, we simply state that our calcula- 
tions include the errors that result from the detector noise, 
sample variance, and the cleaning of CMB using the 150 
GHz channel of EDGE. They do not include the effects of 
contamination by interstellar dust, which we discuss next. 

4.2. Dust 

Galactic dust is far from being an isotropic Gaussian 
random field and is a more dangerous contaminant because 
its frequency dependence is more uncertain and varies spa- 
tially. Moreover, dust can contribute significantly to the 
signal over the entire frequency range where we have a 
hope of measuring the FIRB. A very encouraging point, 
however, is that the amplitude of the dust power spectrum 
varies considerably over the sky. In particular, as we will 
demonstrate, there are a number of sizeable fields with 
substantially less fluctuation power than average, even 
when that average is restricted to high galactic latitude. 
Hence, an important part of any plan to observe FIRB 
fluctuations is to choose fields with very low dust fluctua- 
tion power. 

We begin from the millimeter-wave predictions of the 
two-component dust model of Finkbeiner et al. (1999, 
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Fig. 7. — Histograms of variance in units of the mean variance 
for square patches in an 80° X 80° field centered on the SGP. Light 
dotted lines are for the 256 patches, each 5° X 5°, and the heavy 
solid lines are for the 16 patches, each 20° X 20°. The top panel is for 
dust emission at 300 GHz, and the bottom panel is for synchrotron 
emission at 30 GHz. Dust maps are the predictions at 300 GHz 
(Finkbeiner et al. 1999). The mean dust variances are (9.8 ^K)'^ 
(small patches) and (15 ^K)^ (large patches). Synchrotron maps are 
WOMBAT predictions at 30 GHz. The mean synchrotron variances 
are (22 /iK)^ (small patches) and (40 /iK)^ (large patches). 



hereafter F99). The two-component model has two 
direction-dependent parameters which are fixed by Dif- 
fuse Infrared Background Explorer (DIRBE) and Infrared 
Astronomy Satellite (IRAS) measurements at 100 and 
240 /um. The optical properties of the dust are assumed 
to be spatially uniform and were chosen to give the best 
possible agreement with the COBE/FIRAS data at longer 
wavelengths. Software for making full-sky predictions of 
this model combined with the COBE /DIRBE and IRAS 
data is available from the Wavelength-Oriented Microwave 
Background Anisotropy Team^ (WOMBAT). 

We select an 80° x 80° field centered on the south Galac- 
tic pole (SGP) and calculate the power spectrum at EDGE 
frequency bands in sixteen 20° x 20° sub-fields. All sub- 
fields had power spectra Ci roughly proportional to 1/i^ 
(as found by, e.g., Gautier et al. (1992); Wright (1998)). 
The histogram of the variances of these sub-fields is shown 
in Figure 7. The differences between the fields is very large. 
While the SGP as a whole is separated from the worst of 
the Galaxy's dust emission, 4 of the 16 sub-fields have a 
variance less than 1/16 of the average SGP sub-field vari- 
ance. 

In Figure 6, the "Dust- AVE" lines at each frequency are 
normalized to have a variance equal to the mean variance 
of the SGP field as a whole, whereas the "Dust-LOW" 
curves have a variance 1/16 smaller. The cleaner sub-fields 
have such low dust emission that the FIRB fluctuations 
should completely dominate. In light of this, for the pur- 
poses of forecasting FIRB parameter errors in § 5, we have 
allowed ourselves to treat the dust in a very simple man- 
ner. We assume that the dust emits a grey-body spectrum 
with emissivity index a = 2 and a temperature of 18°K, 
completely coherent between frequency bands, with an an- 
gular power spectrum falling as Ci ex l^^. The amplitude 

** http://astron.berkeley.edu/wombat/ 



is the same as the "Dust-LOW" curves at 320 GHz. This 
single-component model, normalized at 320 GHz, is suffi- 
cient for our purposes in forecasting results from EDGE. 
In the low-dust regions, the dust is sufficiently unimpor- 
tant that our results are insensitive to the details of the 
modeling. To fully understand what can be learned from 
full-sky surveys, such as Planck, we would need a more so- 
phisticated modeling of the dust. Here, we simply assume 
that the "Dust-LOW" amplitudes apply to Planck as well 
but include only 10% of the sky. 

In Figure 6, we also show the results of repeating the 
sub-division of the SGP into 256 sub-fields, each 5° x 5°. 
The smaller patches show even greater variation in the 
dust fluctuation power. Thus, for deep small-field obser- 
vations, one can lower the dust contamination greatly by 
choosing the right fields to observe. This may be particu- 
larly important for attempts to measure the polarization 
of the CMB. 

Finally, for comparison, we have done the same exercise 
with WOMBAT predictions for 30 GHz synchrotron maps. 
These show significant variation, but much less than the 
dust, presumably due to the larger galactic scale-height 
for synchrotron emission. 

5. FORECASTED PARAMETER ERRORS 

With the formalism in place and the physical compo- 
nents specified, we can now proceed to consider how mea- 
surements of the FIRB anisotropics constrain the emissiv- 
ity, bias, and spectral energy distributions (SEDs) of the 
far-infrared sources as a function of redshift. However, we 
will first describe the degeneracy structure of the theory, so 
that the reader can better understand what combinations 
can be constrained and how our quantitative investigation 
should proceed. 

5.1. Degeneracies and Reconstruction 

Two major degeneracies affect the interpretation of 
FIRB anisotropics. First, changes in the emissivity of 
the sources may be compensated by changes in their bias 
(HKOO), as evident in equation (12). Only the product jb 
(eq. [13]) can be constrained by confusion-limited anisot- 
ropy measurements, and so we will work only with this 
quantity. Of course, the emissivity of the sources can be 
measured through the FIRB mean (§2) or through the 
study of sources from higher-resolution observations. Any 
separation of bias and emissivity from anisotropy measure- 
ments will require such external information. 

Second, within confusion-limited observations, changes 
in the bias-weighted emissivity are partially degenerate 
with changes in the spectra of the sources coupled with 
a shift in redshift. Pedagogically, this degeneracy may be 
best understood by first considering the special case of a 
power-law matter power spectrum, for which the degen- 
eracy is complete. 

Imagine that each redshift emits a featureless spectrum 
with a high-frequency cutoff (e.g. a grey-body at some 
temperature). Each redshift will contribute anisotropy 
only to observed frequencies below the redshifted cutoff. 
As one considers higher frequencies, the contributions from 
different redshifts begin to drop out, starting with the 
highest redshifts. By measuring the reduction in anisot- 
ropy as well as the imperfect cross-correlation, one can 
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isolate the contribution of the sources with that apparent 
cutoff frequency. If one knows the rest-frame frequency 
of the cutoff, then one can infer the redshift and use that 
redshift to map the angular power of the emission to its 
spatial power. By dividing the spatial power spectrum of 
the light by that of the mass, one finds the bias-weighted 
emissivity of the sources at each redshift. 

However, if the emitted temperature is unknown, one 
cannot assign a redshift to a particular observed cutoff. 
For example, an increase in temperature of the emission 
would imply a larger redshift. This shift in redshift alters 
the mapping between physical and angular scale; however, 
because the spatial power spectrum is a power-law, the 
shift in distance scale can be compensated by a shift in 
normalization so as to leave the observed angular power 
spectrum unchanged. We thus have a degeneracy between 
the emitted spectrum, the redshift, and the bias-weighted 
emissivity of the sources. In models of interest, a 10% shift 
m 1 + z mimics a 30% shift in fluctuation amplitude (60% 
in power). 

Another way to describe this is to view the FIRB anisot- 
ropics as a sum of correlations of different apparent temper- 
ature components. By analyzing the frequency structure 
of the correlations, one can determine the power spectrum 
associated with each apparent temperature. However, un- 
less one knows the original emitted temperature, one can- 
not determine the redshifts of the apparent temperature 
components. Without the redshift, one cannot convert the 
angular power spectrum to a spatial one. 

The degeneracy is lifted when the power spectrum is 
not a power-law. With a feature or bend in the spatial 
power spectrum, one has a second handle on the redshift 
of the sources through the measurement of where the fea- 
ture appears in angle. By combining the details of the 
frequency covariance with the location of features in the 
angular power spectrum, one can measure the emission 
spectrum separately from the bias and luminosity. Unfor- 
tunately, because CDM power spectra are fairly smooth 
and lack strong features, the redshift-temperature degen- 
eracy is only partially lifted. The slow dependence of the 
comoving angular diameter distance on redshift further 
limits one's ability to measure the redshift. 

The above paragraph assumed that the bias of the 
sources was scale-independent, so that scales in the matter 
power spectrum would shine through in the FIRB anisot- 
ropics. The bias may instead be scale-dependent. Even 
a smooth dependence on scale can shift the location of 
soft bends in the power spectrum, such as are found in 
CDM models. This would confuse the interpretation of 
the angular power spectrum in terms of redshift and skew 
the inferred bias-weighted emissivities. Of course, a gen- 
eral scale-dependent bias will ruin any use of the angular 
power spectrum to measure the redshift of the sources. 

Because of the above concerns, we anticipate that the 
emission temperature and hence redshift of the sources will 
not be well-constrained by the FIRB anisotropics alone. 
In §5.2, we find that biases could be constrained in 4 red- 
shift bins at the level of 5-15% if the emission spectra 
were known, but the large-angle FIRB anisotropics them- 
selves will not deliver the ~ 5% emission temperature con- 
straint that would be needed to break the degeneracy to 
this level. Instead, one must rely on external informa- 



tion, notably from high-resolution infrared imaging with 
foUowup spectroscopy, to supply the emission spectrum as 
a function of redshift. Clearly, the spectral energy distri- 
bution of bright sources will be a good start. We caution 
only that the FIRB anisotropics actually require the bias- 
and luminosity-weighted spectral energy distribution of 
the sources. If the spectra of sources depends strongly 
on luminosity, it may be necessary to extrapolate to un- 
resolved flux levels when integrating over the luminosity 
function. 

Having noted this degeneracy, it is important not to 
over-react to it. First, allowing for uncertainty in the 
emission spectra will make the bias-weighted emission jh 
uncertain, but there will be a combination of these pa- 
rameters that will be well-constrained. Indeed, as shown 
in Appendix A of Eisenstein et al. (1999), the errors on 
the product of jb and a calculable function of the spectral 
variation will be the same as the errors on jb in the limit 
that the spectra are held fixed. In our quantitative work, 
we want to focus on this well-constrained combination of 
parameters; in breaking the degeneracy by fiat, we can 
calculate the errors on this quantity. We will show in §5.2 
that the FIRB anisotropics can deliver several better than 
10% constraints on the properties of FIRB sources. 

Second, the emission temperatures in plausible models 
do not vary by large amounts. If the FIRB anisotropics 
were to reveal a factor of 8 change in emission as a func- 
tion of redshift, it would be unrealistic to explain it by a 
factor of two shift in the temperature scale. While such 
coarse statements "waste" the precision of the FIRB an- 
isotropy measurement, they do show that the measure- 
ments can yield interesting results without detailed exter- 
nal constraints on the emission properties. 

In the same spirit, we will assume that the underly- 
ing cosmology is well known. Introducing uncertainties in 
that sector would cause additional degeneracies. Some are 
simple: the amplitude, shape, or time dependence of the 
matter power spectrum will be directly degenerate with 
the bias. Uncertainties in the relations between distance, 
volume, and redshift will cause more subtle problems akin 
to those described above. However, one can always map 
the constraints on jh into constraints on a product of the 
jb and functions of cosmology. 

5.2. A Quantitative Study 

For a quantitative assessment of what we can learn 
about jb{k,v^z) from measurement of Ci'^ , we split jb 
into four redshift bins with intervals of [0,0.5], [0.5,1], [1,2] 
and [2,4]. We then parameterize jb as follows: 



jb{k, v,z) = ^ Xi{z)atjb{k, t^, z){k/K 



(15) 



where Xiiz) is unity for z in the i**^ redshift bin and zero 
otherwise. The parameters a^ and rii adjust the amplitude 
and scale-dependence, respectively, of the bias-weighted 
emissivity in the i^^ redshift bin. j^{k,iy,z) is simply the 
jb of the HKOO model (and actually has no k dependence) . 
The pivot point fc, is chosen to project to ^ = 500 at all 
redshifts. As we will see, this is nearly the "sweet spot" 
at which the errors on ai and n^ become independent. In 
total, the model has 8 free parameters: 4 amplitudes and 
4 tilts. 
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Table 1 



EDGE Channel Characteristics 
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(^K yi) 
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LI 


150 


2000 


14 


1 




68 


6.2 


79 


< 1 


L2 


220 


1360 


14 


1 




36 


3.3 


76 


< 1 


HI 


320 


940 


6 


6 




91 


9.8 


99 


7 


H2 


390 


770 


6 


6 




71 


7.7 


114 


12 


H3 


480 


630 


6 


6 




60 


6.5 


143 


21 


H4 


590 


510 


6 


6 




49 


5.3 


176 


30 


H5 


710 


420 


6 


6 




42 


4.6 


227 


33 


H6 


870 


340 


6 


6 




36 


3.9 


274 


38 


H7 


1,060 


280 


6 


6 




30 


3.3 


352 


40 


H8 


1,290 


230 


6 


6 




26 


2.8 


454 


39 



NOTES. — Pixel errors apix are for beam-size pixels and observations uniformly covering 
1% of the sky over 10 days. The weight-per-unit solid angle of equation (18) is given by 
w = l/(crpixFWHM^). For the H ("high-frequency") channels, we assume only 4 of the six 
detectors are used. Number of sources observed A''(> 5(7pix) assumes model E (G98). See 
http://topweb.gsfc.nasa.gov for more details. 



We can forecast errors on these parameters by theoret- 
ically propagating the measurement errors. This is conve- 
niently done with the Fisher matrix formalism (Kendall & 
Stuart 1969; Tegmark et al. 1997). In general, the Fisher 
matrix for a set of parameters, Op, depends on the covari- 
ance matrix of the data (in this case the maps) and the 
derivatives of the covariance matrix with respect to the 
parameters: 



J. rt' 



1, 
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_^dC 



c'- 



dC 



da 



pi 



(16) 



where the trace here runs over the suppressed indices £, m 
and u. Calculating the Fisher matrix is most easily done 
in spherical harmonic space where the covariance matrix 
of the maps can be written as 



^ulm^u'l'm' — {P'l^lvaO'u'l' rn' ) — ^l ^ll'^n 



(17) 



with 






Y,Cr'^''^b,{l)K,{l)+w-\,^)5,,, (18) 



where bi,{l) is the window function of the beam, w{i') is 
the weight~per-unit solid angle of the map at frequency v, 
and the index i labels the various components, which we 
have assumed to be uncorrelated. Such multi-frequency, 
multi-component analyses have been performed for CMB 
experiments by a number of authors (Tegmark & Efs- 
tathiou 1996; Bouchet, Gispert & Puget 1996; Knox 1999; 
Tegmark ct al. 2000; Cooray, Hu & Tegmark 2000). We 
use four components in our modeling: the FIRB corre- 
lations, the FIRB shot noise. Milky Way dust, and the 
CMB. We take the shot noise to be uncorrelated between 
different frequencies for the sake of simplicity. In reality, 
the cross-frequency correlations are probably strong. We 
justify our approximation below. 

We can separate the calculation into separate multi- 
poles, so that with uniform coverage over a fraction of 



the sky, /^ky. 
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Tr 
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(19) 



where now all the terms in the trace are matrices with 
(suppressed) indices i^,iy'. When /sky is less than 1, this 
equation is approximate, although it is generally a very 
good approximation. The /sky factor accounts for the de- 
creased number of modes available when the sky coverage 
is less than full. 

At this point, the matrices one needs to invert are only 
size rich by rich) where rich is the number of channels, and 
therefore this calculation can be done quickly. To simplify 
the amplitude derivatives, we take them with respect to 
af instead of a^ . 

The diagonal element Fa gives the variance on the i^^ 
parameter when all the others are being held fixed (the un- 
marginalized case), while the inverse of the Fisher matrix 
yields the predicted covariance matrix when all parame- 
ters are allowed to vary simultaneously (the marginalized 
case). We will also display the results when the amplitudes 
are allowed to vary but the scale-dependence parameters 
are held fixed. For all our calculations we restrict the sum 
in equation (19) to € < 1600 in order to reduce the de- 
pendence of our results on the amplitude of the shot-noise 
power spectra, which are quite uncertain. 

We begin by displaying the predictions for the EDGE 
mission in Figure 8. In the unmarginalizcd case, the two 
parameters of the 3 low-redshift bins can be constrained 
to ^ 1%, while the highest-redshift bin is measured only 
to ^5%. In the marginalized case, all parameters are con- 
strained at the ^ 10% level. Interestingly, the constraints 
on the amplitude are very similar whether or not the tilts 
Tii are allowed to vary. This indicates that the values fc, 
were well-chosen so that the uncertainties in the ampli- 
tudes and tilts are nearly uncorrelated (Eisenstein et al. 
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Fig. 8. — Forecasted relative errors from EDGE observations on 
the ji, amplitudes (bottom panel, squares) and absolute errors on 
the bias spectral index parameters (top panel, triangles) for a model 
with 8 parameters: the amplitude factor a, and n in 4 redshift bins. 
Dotted lines connect results for parameter errors assuming the other 
7 parameter held fixed. Solid lines connect the results for parameter 
errors that follow from assuming none of the parameters are fixed. 
The dashed line shows the result for the amplitude parameters when 
holding only the n parameters fixed. 



1999, for more details, see). 

Since EDGE is covering 1% of the sky to i' « 1000, one 
would have expected it to yield a 1/VO.Ol x 1000^ = 1% 
measurement of some combination of these parameters. 
Indeed, diagonalization of the 4x4 amplitude sub-block 
of the Fisher matrix shows that this is the case. Table 2 
shows that the errors in the fully marginalized case dis- 
played in Figure 8 are dominated by one eigenvalue. Ex- 
cursions in parameter space that move all of the jb am- 
plitudes in the same direction are constrained at better 
than 1%, whereas excursions that move neighboring red- 
shift bins in opposite directions are only constrained to 
10%. Indeed, our choice to use only 4 bins constitutes an 
implicit smoothing prior on the set of all possible excur- 
sions of the function jb{i^, z), and because the FIRB cor- 
relations are projected quantities, one will never constrain 
highly oscillatory star-formation histories. Physically, of 
course, such models are absurd. Models that differ by 
broad, smooth changes can be distinguished by EDGE at 
high accuracy, ~ 1%, which is better than the marginalized 
errors in Figure 8 would indicate. 

We repeat all of the analysis for the Planck mission. 
Since Planck covers the full sky, over much of which the 
dust is a significant contaminant, an accurate treatment 
of how well Planck can measure the a and n parameters 
of equation (15) would require a more careful treatment 
of the dust contamination than has been done here. To 
avoid these complications, we forecast the results for a 
conservative analysis of the Planck data which only uses 
the cleanest 10% of the sky. We assume that the fluctu- 
ation power in this cleanest 10% is the same as assumed 
for the EDGE observations of 1% of the sky. 

The results for Planck are shown in Figure 9 and Table 2. 
The behavior is quite similar. The difference between the 
errors on a^ with and without marginalizing over the Ui are 
larger in the Planck case. This indicates that the ai and 
Ui are somewhat correlated. A slightly larger choice for fc. 
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Fig. 9. — Same as Figure 8 but for Planck. 



would have decorrelated the two sets of parameters, and 
the errors on ai would have been simply those in the case 
where the Ui were not varied. In other words, with a better 
choice of fc* the amplitude constraints would be quoted as 
closer to 5% than 10%. The movement of the sweet spot 
to higher fc is a direct result of the greater sensitivity and 
slightly higher angular resolution of Planck. 

We have studied the dependence of the parameter errors 
to changes in experimental parameters and the parameters 
governing sources of astrophysical noise (namely, the shot 
noise and dust power spectrum amplitudes). Halving the 
dust power spectrum makes less than 10% changes in the 
errors. Increasing it by a factor of 20 makes at most a 
factor of 2 increase in the parameter errors. Increasing 
the shot noise power spectrum by a factor of 2 makes less 
than a 10% change in the errors. 

This robustness is fortunate since, as has already been 
mentioned, our modeling of these contaminants has not 
been very sophisticated. For the dust we have assumed 
complete coherence, and for the shot noise complete de- 
coherence; reality for both is somewhere in between. In 
principle the Op of equation (16) would be not only the 
FIRB parameters, but also dust and shot-noise param- 
eters. In effect, our calculation assumes that the Cj^'^ 
for the dust and shot noise are perfectly known. In real- 
ity these statistical properties will have to be determined 
from the data, although external datasets will provide use- 
ful information as well. A further weakness of our dust 
modeling is the implicit assumption of statistical isotropy. 
A more rigorous treatment of the dust, which one would 
need to study FIRB anisotropy in regions of typical dust 
fluctuation power, would be quite challenging indeed. 

Removal of low frequency channels degrades our high- 2; 
constraints, as we would expect. Eliminating the Planck 
217 GHz channel, while not changing the unmarginal- 
ized errors, does increase the marginalized errors — in the 
lowest-z bin by 50% and the highest-z bin by a factor of 2. 
Removing the three lowest frequency FIRB channels (L2, 
HI and H2) has a similar effect. However, simply remov- 
ing channel L2 makes almost no difference. Presumably 
the importance of LS would be more readily apparent if 
the highest-z bin, which runs from 2; = 2 to z = 4, were 
split up more finely. As the bin stands, the HI and H2 
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Table 2 



DiAGONALIZING THE FiSHER MATRIX 



Experiment 


Error on bias 




Eigenvector 






A-/V2 


< z < 0.5 


0.5 < z < 1 


1< z < 2 


2 < z <4 


EDGE 


0.0040 


0.2799 


0.7094 


0.6451 


0.0471 




0.0151 


0.4525 


0.4860 


-0.7148 


-0.2192 




0.0777 


0.5798 


-0.2458 


-0.0380 


0.7759 




0.1922 


0.6171 


-0.4474 


0.2673 


-0.5896 


Planck 


0.0014 


0.2508 


0.6780 


0.6884 


0.0588 




0.0060 


0.4097 


0.5361 


-0.6469 


-0.3554 




0.0234 


0.4801 


-0.0460 


-0.2023 


0.8523 




0.0912 


0.7340 


-0.5008 


0.2582 


-0.3793 



NOTES. — The eigenvalues and eigenvectors of the Fisher matrices for EDGE and 
Planck. The scale-dependence parameters Ui have been held fixed here, leaving 
only the four amplitude parameters ai. Each row gives an eigenvalue and the 
corresponding eigenvector. The eigenvalues are listed as the inverse square root 
of the Fisher matrix eigenvalue \f divided by two; this is the l-a error on the 
bias. The factor of two comes from our use of af as a parameter instead of a;. 
The eigenvectors specify how that error is divided between the four redshift bins. 
Note that the spectrum of eigenvalues increases rapidly: model changes that would 
move all of the biases up or down together are extremely well constrained (better 
than 1%), whereas changes that shift the biases in alternating fashion are relatively 
poorly constrained (10-20%). 



channels are helping to constrain its amplitude because of 
their sensitivity to just the lower end of the range. 

The high frequency channels are important for sepa- 
rating out the low-redshift contribution from the high- 
redshift contribution. Removing the highest three fre- 
quency channels of EDGE has only a small effect on the un- 
marginalized errors, but the marginalized errors increase 
by factors of ^ 10. For Planck removal of the 857 GHz 
channel results in mostly small changes to the unmarginal- 
ized errors but factors of from 4 to 8 increase in the 
marginalized errors. 

While frequency range is important, we have not been 
able to see, with this study, a benefit from a dense sam- 
pling of that range, such as offered by EDGE. Removing 
all the even EDGE channels (L2, H2, ..., H8) while simul- 
taneously doubling the weight of the odd ones (to keep 
the total weight roughly unchanged) results in practically 
no change in either the marginalized or unmarginalized er- 
rors. This may be due to the fact that our bins of 1 -I- z 
are coarser than the frequency bins; i.e., the finer the fre- 
quency sampling, the faster the changes in emissivity with 
redshift one should be able to detect. Dense sampling also 
makes possible more consistency tests and better monitor- 
ing of possible dust contamination. 

Once again we stress that our results here have assumed 
small departures around a particular model for the FIRB 
fluctuations. In particular, it may be possible to have 
detections of significant emissivity at redshifts beyond 4, 
if there is such emission. Even if there is not such extra 
emission, and Ajh/jh turns out to be large for z > 4, one 
may still be able to set very interesting upper limits on 
jh. Combined with a lower limit on the bias, this could be 



turned into an upper limit on j^ alone — providing a new 
constraint on energy production at high-redshift. 

It is now convenient to give some more details on the 
FIRB power spectrum error forecasts in figure 6. They 
were the same as the parameter error forecasts described 
in this section (i.e, equation (16) through equation (19)) 
with the following changes: (1) for the errors in each panel 
only two channels were used — the one corresponding to 
that panel, and the 150 GHz channel as a CMB monitor, 

(2) the parameters were the FIRB power spectrum multi- 
pole moments themselves, Cj^" , instead of the at and rn, 

(3) dust was ignored, and (4) the FIRB Cj^'^ included con- 
tributions from non-linear corrections. It should be clear 
that the figure 6 error forecasts were only used for that fig- 
ure. They are not an intermediate step in our calculation 
of FIRB parameter errors. 

6. DISCUSSION 

The large-angle anisotropics of the FIRB offer a comple- 
mentary view of the high-redshift universe to that given 
by direct investigation of far-infrared sources. Clearly, 
if one could measure the redshift and spectrum of every 
far-infrared source in some large patch of sky, then one 
would recover all of the information we have described 
and more. However, this will not be possible in the near 
future. By combining high-resolution measurements of 
the average spectrum and emissivity as a function of red- 
shift with the known statistics of large-scale density fluc- 
tuations, the confusion-limited measurement of the FIRB 
anisotropics gives us a window on the large-scale biasing 
of the far-infrared sources. In CDM cosmologies, this bias 
holds implications for the typical mass of the host halo of 
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the sources (Mo & White 1996), thereby giving an obser- 
vational window on the halo environment of the starburst 
and AGN activity that power the FIRB. Theories for how 
the galaxy density-morphology relation arises (KaufFmann 
et al. 1993) and why clusters of galaxies form rather than 
single overmerged galaxies often appeal to there being a 
preferred mass (or velocity) scale for galaxy interactions 
(Kolatt et al. 1999; Somerville et al. 2000). If such inter- 
actions lead to dust-obscured star formation and nuclear 
activity, then the bias of the FIRB sources could inform 
the modeling of bulge and group formation. 

It is interesting to compare the limits on clustering that 
can be gained from confusion-limited measurements to 
those that could be gained from more typical correlation 
analyses of sets of individual sources. 

First, we will consider the case in which exact red- 
shifts for the sources are known and ask what number 
and density of sources would be needed to measure the 
power spectrum over a broad band centered at a comov- 
ing wavenumbcr k = 0.2/iMpc^ {£ « 800) to We assume 
that P{k) = 1000/1^3 Mpc^ at /c = 0.2/iMpc"^; this cor- 
responds to a bias of 2 at z « 2.5. For a fixed number 
of objects TV, the power spectrum error bars are mini- 
mized when the comoving volume density is equal to P^^ 
(Kaiser 1986; Tegmark 1997a); in this case, the fractional 

error is approximately 2v2/a 



TV 



k3p{k) 

271-2 • 



Hence, 5% accu- 



racy requires approximately 8000 sources over a comoving 
volume of 8 X 10^h~^M])c^. ^. If this volume is roughly 
cubic, then this is approximately a range of 0.25 in red- 
shift and 8 square degrees at z = 2.5. If the source density 
is much smaller or larger than 10~^/i^ Mpc~^, then the 
number of sources required to achieve a 5% measurement 
of large-scale power increases dramatically. 

Next, wc consider the case in which the redshifts are only 
approximately known, e.g. by using the far-infrared or ra- 
dio spectral index (CariUi & Yun 1999; Blain 1999e). We 
approximate this situation by saying that the sources can 
be divided into a series of redshift slices. Each slice is an- 
alyzed by angular correlation methods and is presumed to 
be statistically independent from the others. In the limit 
that the conformal distance r changes only slightly across 
the redshift slice, the angular power spectrum C'e of the 
sources in the slice is simply the spatial power spectrum 
at fc = £/r divided by the comoving volume per steradian 
in the slice. For P ~ lOOO/i^"^ Mpc and a slice of unit 
redshift at z = 2.5, this gives Ci ~ 10~^. The optimal 
number of sources per steradian to minimize the errors on 
the measurement of C^ for a given number of sources is 
simply CJ^ . With this surface density, the fractional er- 
rors on d are roughly 2\/2/Jn^^. Inserting I = 800, 

one finds a requirement of roughly 3 x 10"'' sources (in the 
unit redshift band) spread over 100 square degrees in order 
to achieve 5% fractional errors on the power spectrum. 

In both cases, the amplitude of the power spectrum gives 
a preferred surface density in a narrow redshift band. This 
density implies an upper limit to the flux cut of the selec- 
tion. If wc assume that the comoving source density is con- 

Of course, the required number and volume scale as the inverse 
square of the desired accuracy; 1% accuracy would require ^ 200, 000 
sources. 



stant out to z = 5 at this flux level, then integrating over 
redshift imphes ~15,000 total sources per square degree. 
This corresponds to roughly 1 mjy at 850 /im (Hughes 
et al. 1998). One gets the same answer in either of the 
two cases above and regardless of the thickness of the red- 
shift slice in the angular analysis; the value is controlled 
by the amount of volume at high-redshift and the ampli- 
tude of the power spectrum one is trying to measure. The 
detection of sources becomes confusion-limited when the 
inverse of the source surface density is less than 30 times 
the beam area. Hence, one can immediately calculate that 
detecting sources at the density needed to maximize one's 
measurement of the high-redshift power spectrum with a 
fixed number of sources requires a beam smaller than 5" . 
Surveys with beams larger than this will reach their con- 
fusion limit before achieving the optimal source density 
and will need to measure more sources in order to gain 
the same fractional limits on the power spectrum. Sur- 
veys with beams smaller than this should not integrate to 
the confusion-limit but instead cover more sky. 

Clearly, resolving the FIRB sources has other benefits 
with respect to clustering that cannot be matched by 
a confusion-limited measurement; for example, one can 
study the clustering of different populations of sources. 
However, the above calculations show that achieving ac- 
curacy equivalent to Figure 8 on the large-scale power 
spectrum requires large sets of sources: hundreds of thou- 
sands with crude redshifts or many thousand with accurate 
redshifts. 

An aspect of the FIRB anisotropics that we have not 
included in our treatment is the possibility of stochastic 
bias (Dekel & Lahav 1999). Our assumption that the bias 
between the emission at different frequencies is perfectly 
correlated (i.e. that the frequency covariances can be de- 
composed into b{v)b{v')) is optimistic because it allows 
measurements of the cross-correlations between frequen- 
cies to constrain the jh directly. However, we find that 
the HKOO predictions of the cross-correlation differ suf- 
ficiently little from perfect correlations that the leverage 
on the redshift decomposition is relatively weak. Break- 
ing the assumption would slightly weaken our constraints 
but would also open up a new sector of physical parame- 
ters, namely the stochasticity of the relative bias between 
different frequencies. This can occur if at a given red- 
shift, there are multiple components with differing spectral 
energy distributions that are imperfectly correlated with 
each other (and thus necessarily with the density). Such 
flawed correlations are necessary in the non-linear regime, 
where the density contrast can't fall below —1, but it is un- 
known whether bias is stochastic on the large scales stud- 
ied here. If the observed correlation coefficients were more 
distinct from unity than the baseline predictions, the lack 
of correlation could be measured and perhaps interpreted. 
Obviously, it would be intriguing to associate the sepa- 
rate populations with active nuclei versus star formation 
or quiescent star formation versus starbursts. Hence, the 
cross-correlations between frequencies could be an interest- 
ing route for FIRB anisotropics to constrain the properties 
of high-redshift sources. 

There remains the possibility of truly diffuse emission 
contributing to the FIRB, i.e., emission which does not 
resolve into distinct sources. As an example, grey dust 



14 



distributed throughout the IGM has been invoked as a 
non-A explanation for the anomalously dim high- 2; super- 
novae (Aguirre 1999), and it has been shown that if this 
explanation is correct, then a significant contribution to 
the FIRB at low frequencies comes from this grey dust 
(Aguirre & Haiman 2000). We can expect the cluster- 
ing properties of such emission to be different from that 
of the individual sources. A cross-correlation analysis of a 
wide-field map of resolved FIRB sources with a map of the 
unresolved FIRB emission (as measured by, e.g., Planck or 
EDGE) could eventually be used as a probe of any such 
diffuse component. 

7. CONCLUSIONS 

We have presented a study of the large-angle cluster- 
ing properties of the far-infrared background, focusing on 
the role of confusion-limited anisotropy experiments. We 
find that under reasonable assumptions, and in regions of 
the sky with unusually low Galactic dust emission, cor- 
related FIRB fluctuations should dominate the observed 
angular power spectrum over a wide range of frequencies 
and angular scales. As a result, such experiments will pro- 
vide useful constraints on any model for the sources of the 
FIRB. The main drawback of the confusion-limited tech- 
nique is that it cannot easily separate different physical 
classes of sources, and so many important questions about 
the nature and distribution of far-infrared sources will re- 
main squarely in the domain of high-resolution studies. 
However, matching the constraints on the large-scale clus- 
tering that could be produced by the EDGE or Planck ex- 
periments requires an impressively large set of point source 
measurements with or without accurate redshift informa- 
tion. 

To model the FIRB fluctuations, one must specify the 
spectrum, emissivity, and bias of high-redshift galaxies 
and AGN, as well as the redshift evolution of these quan- 
tities. Despite the inherent degeneracies, we find that 
confusion-limited anisotropy observations yield tight con- 
straints on certain combinations of these model ingredi- 
ents. Without external spectral information, the emis- 
sivity is degenerate with the temperature of the sources. 
However, as the plausible ranges of temperatures and bi- 
ases are reasonably small, the FIRB anisotropics by them- 
selves could detect strong variations in the emissivity as 
a function of redshift, notably the turn-on of embedded 
star formation at high redshift. We find that if the mean 
emission spectrum was precisely known, the bias-weighted 
emissivity could be measured to ~ 10%; while strong priors 
on the amount of oscillation allowed as a function of red- 
shift can allow these constraints to approach 1%. These 
constraints are tighter than differences between possible 
models for the sources of the FIRB, rendering such stud- 
ies of the infrared sky a useful probe of the high-redshift 
universe. 
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